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New error bounds are given for approximate solutions of differential equations of the form 

d'wId'C - {u 2 (±a 2 -n + f(u 9 a,$)}w 

in terms of parabolic cylinder functions. Here f is a real variable, a is a bounded real parameter, 
a is a large positive parameter, and \p(ii,a^) is a continuous function of a and ?. Zero values of 
a are admitted. 
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1 . Introduction and Summary 

Recently [2] 1 , [3], I developed a uniform asymptotic theory of second-order linear differential 
equations with two coalescing simple turning points, and applied the results to the associated Legendre 
equation. Subsequently, in the course of writing a paper on connection formulas for multiple turning 
points [4] it became clear how to effect some improvements in two of the four cases treated in [2] . The 
purpose of the present note is to describe these improvements. It will be assumed that the reader is 
familiar with the results presented in [2], and the same notation will be used except where indicated 
otherwise. 

The next section introduces new auxiliary functions for the solutions of the modified Weber equa- 
tion. The new form of the general approximation theorem is stated and discussed in the third (and con- 
cluding) section. 

An application of the results to the approximation of Whittaker functions with both parameters 
large will be published in due course. 



2. Auxiliary Functions for the Modified Weber Functions 

The modified Weber functions W(b,x) and W(b, — x) are solutions of the differential equation 

d 2 w/dx 2 = {b — \x 2 )w. 
As in [2] , §8.6 we denote by o-(6) the smallest positive root of the equation 

k~ 1/2 (b)W(b,x) =k 1/2 (b)W(b,-x), 



* Present address : Institute for Physical Science and Technology, University of Maryland, College Park, Md. 20742. 
1 Figures in brackets indicate the literature references at the end of this paper. 
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but we now extend this definition to all real values of the parameter b. We already know that 

(7(6) ~26 1/2 , 6->+ 00. 

Next, because <r(b) is bounded by the smallest positive zero of W(b,x), the expansion (15.7) of [1] 
shows that 2 

a(b) -»0, &-»— 00. 

In place of the definition [2] , (8.15) of the weight function E(b,x) , we now define 

/ k 1/2 (b), - 00 <x< -*(b), 
£(&,%)=: 1 {W(b-x)/W(b,x)} 1/2 ,-<r(b) <x<cr(b), 
( k~ 1/2 (b), <j(b) <X< 00. 

Then E(b,x) is continuous for all real values of b and x, and is also a nondecreasing function of x. 
Furthermore, denoting 1/E(b,x) by E" 1 (b,x) , we have 

E[b—x) —E-^ib^x). 

Using the above definition of E(b,x) we now define modulus and phase functions for all real values 
of b and x by the relations 

W(b,x) — E-Hb^Mikx) sinO{b,x), 

W(b—x) — E(b,x)M{b,x) cos 6(b,x), 

W'{b,x) — E^ib^Nib.x) $mw{b,x), 

W'(b—x) = E(b,x)N(b,x) cosa>(6,%). 

Thus when — co < x < — o-(b) we have 

M(b,x) = {k(b)W 2 (b,x) +k-Hb)W 2 (b-x)} 1/2 , 

0(b,x) = tan-^k^Wib^/Wib-x)}, 

N(b,x)= {k(b)W' 2 (b : x) +k- 1 (b)W' 2 (b-x)} 1/2 , 

<o(b,x) = tan-^kWW'ib^/W'ib-x)}. 

Secondly, when — o-(6) <#< a{b) we have 

M(b,x) = {2W(b,x)W(b-x)} 1/2 , 0{b,x)-l^ 

N(h , _ \ W 2 (b,x)W 2 {b,-x) +W' 2 (b-x)W 2 (b,x) \v 2 . 
iv (a*) - <^ W(b,x)W(b-x) ) 



(h s t _J W'(b,x)W(b-x) { m 
„(&,*) = tan il^^—^^-V 



Lastly, when <r{b) < x < oo we have 

M{b,x) = {k-Hb)W 2 (b,x) +k(b)W 2 (b,-x)} 1 <\ 
6{b,x) = lan- l {k- l {b)W{b,x)/W{b-x)}, 
N(b,x)= {k- l {b)W' 2 {b,x) + k(b)W' 2 (b-x)} 1/2 , 
„,(&,*) = tan-^k-H^W'^x) /W'{b,-x)}. 



2 The actual asymptotic form of cr(6) as b -» — oo can be shown to be &(b) , — ' ^|fe| -1 / 2 e x& , but we shall not make use of this more precise 
result in the present paper. 
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Each of the functions M(b,x), 6(b,x), N(b,x), and w(b,x) is continuous in b and x. Reflection 
properties are given by 

M(b-x) = M(b,x), 0(b—x) =^ — 0(6,*), 

N(b—x) = N(b,x), uib—x) = \ir — co(b,x) . 



As x — » ± oo 



M(b,x) ~ \2/x\ 1/2 , N(b,x) ~ \x/2\ 1/2 . 



Compared with [2], §§8.6 and 8.7, the modulus functions M(b,x) and N(b,x) agree with their 
previous definitions, except when b < and — o-(b) < x < o-(6). On the other hand the phase func- 
tions agree only in the case of 6(b,x) when b > and x > — a(b) . 

3. Approximation Theorem 

The following theorem replaces both Theorems III and IV of [2] . The definitions of E, M, and /V 
are those just given in §2. 

Theorem: Assume that for a e [0,A], f € (£i,fo), and each value of the positive parameter u, the func- 
tion i/f(u,a,£) is continuous and <7/u,£2(F) is bounded. Then in this region the differential equation 

d 2 w/d£ 2 = (u 2 (=ra 2 -£ 2 ) + ^(u,a,£) }w (3.1) 

has a continuous solution w(u,a,£), with continuous first and second partial £ derivatives, given by 

w(u,a,£) = W(±iua 2 ,£V2^) + e(u,a,£), (3.2) 

where 

|i(u,a,t)|_ f |3e(u,a,£)M| _ 

M ( ±iua 2 A/2u) (2u) 1/2 N ( ±iua 2 ,fV2u) 

< E-'(±iua 2 ,gV2u)[exp| n *ffiJ Q/ t , fa (F)|- l], (3.3) 



and 



1(b) = sup {n(x)M 2 (b,x). } (3.4) 

x e( — oo, oo ) 



In this result either the upper sign or the lower sign is to be taken throughout. Upper signs cor- 
respond to Theorem III of [2], and lower signs correspond to Theorem IV of [2]. In hoth cases a 
numerically satisfactory companion solution to w(u,a£) may be obtained by applying the theorem to 
the equation obtained from (3.1) by replacing £ by — £. The proof of the theorem is similar to the 
proofs of Theorems I, II, III, and IV of [2], and it is unnecessary to record details. 

The new result simplifies the general theory by replacing two theorems by a single theorem. The 
main advantage, however, is that error bounds are now supplied directly for all real values of f. In 
contrast, Theorem III furnishes error bounds only for nonnegative values of £; compare [2], §9.4. 
Theorem IV, however, supplies error bounds for all real values of £, and a property of Theorem IV that 
is not shared by the new theorem is that it is possible to prescribe anywhere in the closure of (£1,^2) a 
reference point £0 at which e(u.a^) and de(u.a^) /df both vanish. In consequence, it is conceivable that 
there may be applications in which Theorem IV is to be preferred. The new theorem should generally 
be used in place of Theorem III, however. 

It should be noted that although the structure of the error bounds (3.3) is independent of the 
choice of ambiguous sign in (3.1), (3.2), and (3.3), the asymptotic nature of the bounds as u -> co 
differs in the two cases. This is in consequence of the results 

1(b) =0(6 1 /3) j 6->+ 00, (3.5) 

1(b) =0(1), 6->-oo, (3.6) 
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which correspond to similar results in §§9.2 and 10.2 of [2], Thus, for example, if |^(w,a,£)/f| is inde- 
pendent of u and integrable at f = — co and -f- co uniformly with respect to a, then by taking the 
balancing function Cl(x) in the definition [2], (6.3) of the error-control function F{u,a£) to be 
1 -f- \x\, we may show that 

e(u,a,Z) = E-Hiua^^~2^)M(^ua 2 ^2ii)0(u- 2/ nn u), (3.7) 

when upper signs are taken, compared with 

e(u,a,£) — E-H—^Ua^^^Mi-iua^tV^) ^- 1 !!!^, (3.8) 

when lower signs are used. The relations (3.7) and (3.8) apply as u — » oo and are uniform with respect 
to a e [0,A] and f e ( — oo,oo ) ; compare [2], §§9.3 and 10.3. Similarly for the derivative de(u,a,f)/df. 



This research has been supported in part by the U.S. Army Research Office, Durham, Contract 
DA ARO D 31 124 73 G204, and the National Science Foundation, Grant GP 32841X2. 
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